Optical absorption and activated transport in polaronic systems 
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We present exact results for the optical response in the one-dimensional Holstein model. In 
particular, by means of a refined kernel polynomial method, we calculate the ac and dc electrical 
conductivities at finite temperatures for a wide parameter range of electron phonon interaction. 
We analyze the deviations from the results of standard small polaron theory in the intermediate 
coupling regime and discuss non-adiabaticity effects in detail. 
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I. INTRODUCTION 

The investigation of transport properties has been 
playing a central role in condensed matter physics for 
a long time. In recent years optical spectroscopy, for 
example, has contributed a lot to unravel the complex 
physics of highly correlated many-body systems, such 
as the one-dimensional (Id) MX chainsji the quasi-2d 
high-temperature superconducting cuprates^ or the 3d 
colossal magneto-resistive manganites.- That way, op- 
tical measurements proved the importance of electron- 
phonon (EP) interactions in all these materials and, in 
particular, corroborated polaronic scenarios for model- 
ing their electronic transport properties at least at high 
temperatures 4^ 

Polarons are quasi-particles composed of an electron 
and the surrounding ions which in a polar solid, pro- 
vided the electron lattice interaction is sufficiently strong, 
are displaced from their equilibrium positions due to the 
presence of the electron. This bootstrap relation be- 
tween electron and lattice displacement makes the parti- 
cle heavy, because it has to drag with it the potential well 
of the phonons. Polaron motion is largely understood and 
has been worked out theoretically in two important lim- 
its: In the first case the electronic bandwidth is large and 
there is a only a slight change in the particle's effective 
mass due to the EP coupling. These quasi-particles are 
called large polarons or Frdhlich polarons. In the second 
case it is assumed that the bandwidth is small, whereas 
the EP interaction is strong and short-ranged. Now pola- 
ronic effects trap the electron at a certain lattice site and 
the size of the quasi-particle becomes comparable to the 
inter-atomic lattice spacing. Thermally activated hop- 
ping will necessarily be the dominant transport process of 
such small or Holstein-type polarons. Although there are 
experimental systems with clear large and small polaron 
characteristics, most of the above-mentioned novel ma- 
terials belong to the transition region between these two 
limiting cases. Here the relevant energy scales are not 
well-separated and perturbative approaches cannot de- 
scribe the complicated transport mechanisms adequately. 

A recent non-perturbative dynamical mean-field study 



of the Holstein model in infinite dimensions 7 reports 
quantitative discrepancies of the temperature depen- 
dence of the resistivity from standard polaron formulas, 
but the dynamical mean-field approach inherently does 
not account for vertex corrections to the conductivity and 
for longer ranged hopping processes induced by the EP 
interaction. On the other hand, exact numerical investi- 
gations of the Holstein model provided quite a number of 
reliable results for the zero-temperature optical conduc- 
tivity in the Id and 2d (extended) Holstein models. 8 ' 9 

Motivated by this situation, in the present work we 
use a recently developed extension of the Kernel Poly- 
nomial Method (KPM)^i2*ii a refined numerical Cheby- 
shev expansion technique, to compute both the dc and 
ac hopping conductivities at finite temperatures without 
any serious approximation. 



II. MODEL AND METHOD 

Our starting point is the Id tight-binding Holstein 
HamiltonianjiS 

H = - t J2 c h ~ \/^I>* + b i K + fo"o b t b i ' 

(1) 

describing a single electron coupled locally to a 
dispersion-less optical phonon mode, where c\ (b\) de- 
notes the corresponding fermionic (bosonic) creation op- 
erator, and rii = c' i c i . Setting fi — 1 and measuring all 
energies in units of the nearest-neighbor hopping integral 
t, the physics of the model is determined by the dimen- 
sionless EP coupling constants 



A = e p /2t and g 2 — e p /ujq 



(2) 



in the adiabatic {u>o/t < 1) and anti-adiabatic (u>o/t 3> 
1) regimes, respectively. In Id the crossover from large 
to small polaron behavior takes place at A ~ 1 (g 2 ~ 1) 
in the former (latter) case»iii4 Whether large polarons 
form in the Holstein model for d > 1 is still under 
debate, 15 i 16 i 17 
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Addressing the linear response of our system to an ex- 
ternal (longitudinal) electric field we consider the Kubo 
formula for the electrical conductivity at finite tempera- 
tures, which isi& 



Re ct{uj) = 7r 



m.n 



-f)E n _ & -/3E n 



ZLuj 



\(n\j\m)\ 2 8(w - w„ 



,)• 

(3) 

Here Z — e P En is the partition function and 

(3 = T _1 denotes the inverse temperature. Since the 
Holstein Hamiltonian involves bosonic degrees of free- 
dom, the Hilbert space even of a finite L-site system 
has infinite dimension. In practice, however, the con- 
tribution of highly excited phonon states is negligible at 
the relevant temperatures, and the system is well ap- 
proximated by a truncated phonon space with at most 
M(A, g, ojq; T) phononsiiS Then |n) and \m) are the eigen- 
states of H within our truncated D-dimensional Hilbert 
space, E n and E m are the corresponding eigenvalues, and 
Lu mn = E m — E n . In Eq. @ the current operator has the 

standard hopping form, j = i e£^i( c l c i+i ~ c l+i c i) > an d 
connects states with different parity^ Thus, assuming 
that the ground state is non-degenerate, the expectation 
value (0|J|0) vanishes in the absence of an external elec- 
tric field. The limit uo — > of JSJl yields the dc conduc- 
tivity, whereas the optical absorption is given by the fi- 
nite frequency data. For small polarons both results are 
interesting and have been evaluated analytically at an 
early staged At T = 0, the regular part of the optical 
conductivity, cr rog (w) = £ J2n>o w n0 \( n \j\°) \ 25 i u ~ w «o), 
was calculated for finite Id and 2d lattices with periodic 
boundary conditions (PBC) in a wide parameter range of 
the Holstein model, using a combination of the Lanczos 
algorithm and the KPM. 8 

At finite temperatures, a similar straight-forward ex- 
pansion of the conductivity is spoiled by the presence of 
the Boltzmann factors and the contribution of all matrix 
elements between eigenstates of the system. Instead, it 
turns out that a new generalised KPM schem o 11 ' 22 can 
be based upon a current operator density 
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\(n\j\m)\ 2 S(x-E n ) S(y - E„ 



(4) 



Being a function of two variables, j(x, y) can be expanded 
by a two-dimensional KPM, 
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y^ 1 mwki9k9iTk(x)Ti(y) 
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(5) 



where the tilde refers to a rescaling of energy (H — > H ) 
that maps the spectrum of H into the domain [— 1, 1] 
of the Chebyshev polynomials of first kind Tj~(x). The 
finite order N of the expansion leads to Gibbs oscilla- 
tions which can be damped by introducing appropriate 
damping factors. Here we use derived from the Jack- 
son kerneiifli 2 ^, and N = 512 throughout the paper. 
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FIG. 1: (Color online) Renormalized current operator density 
j(x, y) used in the 2d KPM. Data obtained for the Id Holstein 
model with A = 0.2 (left panel) and A = 2.0 (right panel) at 
uio/t = 0.4 (L = 6, M = 50; PBC). 



The core of the numerical work is the iterative calcu- 
lation of the moments fi^i = Tr(Tfe(i/)jT;(_ff)j), where 
the trace can be replaced by an average over a relatively 



small number of random vectors 



Finally, the factors 



l/ w fei = (2 — 5fco) (2 — <5zo) account for the correct normal- 
isation. Given the operator density j(x,y) we find the 
optical conductivity by integration 



Rea(uj) = 



ZLuj 



-/3y. 



dy. 



(6) 

The partition function Z = J_ p(E) exp(—(3E) is eas- 
ily obtained by integrating over the density of states 
p{E) = J2n=o ^(-^ ~ E n ), which can be expanded in par- 
allel to j(x, y). Note the main advantage of this approach: 
The current operator density that enters the conductivity 
is the same for all temperatures, i.e., it needs to be ex- 
panded only once. Figure ^ exemplifies the different form 
of j(x,y) in the weak and strong EP coupling regimes. 

At very low temperatures, the numerical evaluation 
of expression 10 requires some caution, since the Boltz- 
mann factors heavily amplify small numerical errors in 
j(y + uj,y). We can avoid these problems, occurring 
mainly at the lower bound of the spectrum, by treat- 
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FIG. 2: (Color online) Schematic setup for the calculation 
of the finite-temperature optical conductivity (left panel). 
Lowest eigenvalues of the Holstein Hamiltonian for L = 6, 
M — 25, and PBC (right panel). The shaded area marks the 
six lowest eigenvalues to be separated from the rest of the 
spectrum. 
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FIG. 3: (Color online) Optical conductivity in the Id Holstein 
model at T = (in units of ne 2 t 2 ) compared to the analytical 
small polaron result Eq. ® [dashed blue lines]. Exact diag- 
onalization data (ED) are obtained for a system with L — 6 
and M = 45; (Jo is determined to give the same integrated 
spectral weight of the lu > (regular part) of Recr. 



FIG. 4: (Color online) Optical absorption by Holstein po- 
larons at finite temperatures in the adiabatic regime (L = 6, 
M — 45). Dashed curves give the analytical result for finite- 
temperature small polaron transport i-a— The deviations ob- 
served for high excitation energies at very large temperatures 
are caused by the necessary truncation of the phonon Hilbert 
space in ED. 



ing the contributions of the ground state and. some of 
the lowest excitations separately. This is illustrated in 
Fig.0 Wc split the optical conductivity into three parts, 



Rea(w) =RecT tu (w) + Rea^p M (w) + ReCT^ PM (w), (7) 

where the first contribution describes the transitions (ma- 
trix elements) between the S separated eigenstates, the 
second part those between the separated states and the 
rest of the spectrum, which can be expressed as stan- 
dard Id KPM expansions, and finally the transitions 
within the remaining D — S states of the spectrum, han- 
dled by a 2d expansion. Using the projection operator 
P = 1 — y^ a _Q |s)(s|, the moments for the contributions 
RecTi d (w) read ^ = {n\jPT k (H)Pj\n). Of course, the 
number of states one has to separate depends on the 
physical situation. The right panel of Fig. [21 gives the 
lowest eigenvalues of the Holstein model at various cou- 
pling strengths. In the strong-coupling regime (A ^> 1) 
states belonging to the lowest small polaron band have 
almost the same energy as the ground state and therefore 
should be treated separately (cf. the curves for A = 1.5 
and 2). Obviously the situation is far less dramatic at 
weak EP couplings. 
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III. NUMERICAL RESULTS AND DISCUSSION 

A. ac concductivity 

We now apply our numerical scheme to the calculation 
of the optical absorption in the Id Holstein model. The 
results for Recr(a;) and possible deviations from estab- 
lished polaron theory are important for relating theory 
and experiment. The standard description of small po- 
laron transport^ 2 - 5 provides the ac conductivity at T = 



Reo-(w) 



Co 



■ exp 



{u - 2e P ) 2 

4£T p W 



(8) 



For sufficiently strong coupling this formula predicts a 
weakly asymmetric Gaussian absorption peak centered 
at lo = 2s p . A similiar analytical formula can be derived 
for finite temperaturesii&Si- 

Starting at zero temperature, Fig. [3] shows Re<j(a;) for 
various EP-parameters. For A = 2 and coo/t = 0.4, i.e., 
at rather large EP coupling, but not in the extreme small 
polaron limit, we find also a pronounced maximum in the 
low-temperature optical response, which, however, is lo- 
cated, somewhat below 2e p = 2g 2 ui , being the value for 
small polarons at T = 0. At the same time, the line- 
shape is more asymmetric than in standard polaron the- 
ory, with a weaker decay at the high-energy side, which 
fits even better the experimental behavior observed in po- 
laronic materials such as Ti02^ Varying the parameters 
significant discrepancies to a Gaussian-like absorption are 
found. Then the polaron motion is not adequately de- 
scribed as hopping of a self-trapped carrier almost local- 
ized on a single lattice site. 

At finite temperature two different transport mecha- 
nism can be distinguished. Clearly, coherent transport, 
which for large EP couplings is related to diagonal (zero- 
phonon) transitions within the lowest extremely narrow 
polaron band, will be negligible at high temperatures. 
For instance, the amplitude of the current matrix el- 
ements between the degenerate states with momentum 
K = ±7r/3 (A" = 0, ±7r/3, ±27r/3, and ir are the allowed 
wave numbers of a 6-site system with PBC) is of the or- 
der of 10 only. Whereas phase coherence is maintained 
during a diagonal transition, the particle loses its phase 
coherence if its motion is triggered by (multi-) phonon 
absorption and emission processes. These so-called non- 
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diagonal transitions which, of course, can take place also 
with zero energy transfer, become more and more im- 
portant as the temperature increases. Accordingly the 
main transport mechanism is thermally activated hop- 
ping, where each hop becomes a statistically independent 
event. In the small polaron limit, where the polaronic 
sub-bands are roughly separated by the bare phonon fre- 
quency (cf. Fig. |21 right panel), this happens for T > ujq. 
Let us consider the activated regime in more detail (cf. 
Fig. 21 With increasing temperatures we observe a sub- 
stantial spectral weight transfer to lower frequencies, and 
an increase of the zero-energy transition probability in 
accordance with previous results^ 

In addition, we find a strong resonance in the ab- 
sorption spectra at about to ~ 2t, which can be eas- 
ily understood using a configurational coordinate pic- 
ture. Placing a homogeneous lattice distortion u at 
L u consecutive sites by applying the unitary trans- 
formation S^(u) = Ylf " Sj(u) = J|f " exp[u(&| — 
b i )] 1 the transformed Holstein Hamiltonian takes the 
form if = (0\SHu)HS(u)\0) ph = -i£ (i>j) c\ Cj - 

2^/E p oJou n i + <j->ou 2 L u . In the adiabatic strong-EP- 
coupling regime, the ground state can be approximated 
as an electron localized at a certain single site with 
the lattice being in a shifted oscillator state (L u = 1). 
That is, |* ) = |l)eJ O S\(g)\0) ph and E = -e p , in 
accordance with the Lang-Firsov approximation. Now 
let us consider excitations from this ground-state, where 
the lattice distortion spreads over two neighboring sites 
(L u = 2) and the electron is in a symmetric or an- 
tisymmetric linear combination of |l) e ; and |2) e ;, i.e., 
the particle is mainly located at sites 1 and 2 but de- 
localized between these sites. We then find |\&i,±) = 
(\l) el ±\2) el )®St(g/2)Sl(g/2)\0) ph and£i >± = T*-e p /2. 
Whereas the (potential) energy related to the displace- 
ment field is reduced to — e p j L u , the kinetic energy comes 
into play since hopping processes between 1 and 2 are 
allowed. The current operator j connects these different- 
parity states with perfect overlap |(\f r i )+ | - j|'I'i ) _)| = (et) 2 , 
giving rise to a strong signal in the optical absorption. 
Note that the excitation energy = 2t is indepen- 

dent of e p . In order to activate these transitions ther- 
mally, the electron has to overcome the "adiabatic" bar- 
rier A = Ei + — Eq = e p /2 — t. A finite phonon frequency 
will relax this condition. From Fig. we find the signa- 
ture to occur above T > 0.5i. Obviously this feature is 
absent in the standard small-polaron transport descrip- 
tion which essentially treats the polaron as a quasiparti- 
cle without resolving its internal structure. Owing to the 
infinite number of neighboring sites it is also absent in 
the DMFT calculation. Of course, one could also extend 
this scenario to excitations where the electron is delocal- 
ized over more than two distorted lattice sites, but for the 
present parameters the signature of these weakly bound 
states would be rather small. 

Entering the non-adiabatic regime of large phonon fre- 
quencies at fixed A = 2, the pattern of sub-bands sepa- 
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FIG. 5: (Color online) Optical absorption by Holstein po- 
larons at finite temperatures in the non-adiabatic regime 
(L = 6, M = 30). Note that now the abscissa is scaled with 
respect to the phonon frequency. 



rated roughly by uj becomes more pronounced, but is, 
of course, washed out at higher temperatures (see Fig. [S] 
for u>o/t — 0.8). In Fig. the average number of phonons 
contained in the ground state (oc g 2 ) is smaller (g 2 = 5) 
than in the previous case where g 2 = 10 (u>o/t = 0.4). 
This also concerns the activated region ujo/t < T < A 
but for these parameters the simple adiabatic picture an- 
ticipated above breaks down anyway. 

Now let us decrease the EP coupling strength at small 
phonon frequencies u>o/t — 0.2 keeping g 2 = 10 fixed. 
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FIG. 6: (Color online) Optical absorption in the adiabatic 
intermediate EP-coupling regime (the notation is the same as 
in Fig. |1] again we use L = 6, M = 40). The inset illustrates 
the finite-size dependence of Recr(cj) for T ~ (the vertical 
dotted line gives the phonon absorption threshold) . It demon- 
strates that the gap observed at low frequencies and tempera- 
tures is clearly a finite-size effect, i.e., at weak-to-intermediate 
couplings the discrete electronic levels of our finite system 
show up in the conductivity spectra. These effects, of course, 
are of minor importance at larger EP couplings, where the 
polaronic bandwidth is strongly reduced, as well as for high 
temperatures. 
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Results for the optical response in the vicinity of the 
large to small polaron crossover (A = 1) are depicted 
in Fig. Here the small polaron maximum has almost 
disappeared and the 2t-absorption feature can be acti- 
vated at very low temperatures (A — > for the two-site 
model with A = 1). The overall behavior of Rea(u) re- 
sembles that of polarons of intermediate size. At high 
temperatures these polarons will dissociate readily and 
the transport properties are equivalent to those of elec- 
trons scattered by thermal phonons. Let us emphasize 
that many-polaron effects become increasingly important 
in the large-to-small polaron transition region. 29 As a re- 
sult, polaron transport might be changed entirely com- 
pared to the one-particle picture discussed so far. 



B. Sum rules 

Before we consider the temperature dependence of the 
dc conductivity, it is useful to test the sum rules for the 
real part of the optical response. First we have the so- 
called /-sum rule, 



/ Rea(u>)du> = — Eki n ■ 

J-oo L 



(9) 



which relates the w-integrated Re cr(u>) to the kinetic en- 
ergy E kin = -tY,{u 3 )( c l c j)T- Note that the u) = 
(Drude) contribution is included in (0 . The second sum- 
rule for Refi(w) is 
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(jj'Rea(uj)duj = —{j 2 )t- 



(10) 



FIG. 7: (Color online) Conformance of the /-sum rule © in 
the weak (upper panel) and strong (lower panel) EP coupling 
regimes of the Id Holstein model. Results for Stot, partial 
fitot's an d -Ekin are given by solid, dot-dashed and dashed 
lines, respectively. For further explanation see text. 



Throughout our calculations Eq. was fulfilled within 
numerical accuracy, where the thermal average (j 2 )t was 
determined again using a Id KPM. Testing sum rule © 
we gain important information about finite-size effects. 

Figure [7] compares £"kin and StotL/ne 2 obtained from 
our finite-cluster calculation. At weak EP couplings the 
kinetic energy is a strictly monotonic increasing func- 
tion of temperature and becomes strongly suppressed 
at high temperatures due to scattering of the electron 
by thermal phonons. Whereas the /-sum rule is almost 
perfectly fulfilled for smaller values of (3, we found pro- 
nounced deviations at low temperatures which, without 
doubt, can be assigned to the finite size of our Holstein 
ring (cf. the L dependence of S to t shown in the upper 
panel of Fig. and also the discussion of the dc con- 
ductivity below). Clearly finite-size effects become im- 
portant when the temperature is comparable to the en- 
ergy gaps in the spectrum of H (see Fig. EJ. At strong 
EP couplings the transport is hopping-dominated and 
the kinetic energy exhibits a maximum at a finite tem- 
perature that can be related to the thermal activation 
energy of polarons. Now, in view of the narrow small 
polaron band (cf. Fig.0), finite size gaps are small and 
the /-sum rule is fulfilled down to very low temperatures. 



Small deviations appear to vanish rapidly with increas- 
ing system size (see inset, lower panel of Fig. 0). In or- 
der to analyze the contributions of different transport 
processes to -Ekin in some more detail, we have decom- 

^S-l c ld(s) 



posed S tot = S£? + EsJo #tot W + Stot in analogy to 
Eq. 0. Figure [3 shows that the coherent (intra-band) 
contribution (cx «S™) is almost negligible at the tem- 
peratures considered. Inter-band transitions connecting 
eigenstates of the lowest polaron band to higher excited 
states (Id KPM) are the determining factor at low T. 
If the renormalised polaron bandwidth is small enough, 
all states in the band are equally populated, leading to 
pretty much the same values of • At high tempera- 
tures, of course, the transitions covered by the 2d KPM 
are predominant. 



C. dc conductivity 



The dc conductivity is obtained by taking the limit 
-> of ©, which with lim^ (l - e~^)/uj = (3 yields 
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FIG. 8: (Color online) Top row: DC conductivity as a func- 
tion of temperature at weak (left panel) and strong (right 
panel) EP coupling. The bottom panel shows, for strong cou- 
pling, a comparison with recent DMFT results- and standard 
polaron theory-^. 



polaron hopping transport manifests in the temperature 
dependence, leading to the well-known absorption max- 
imum in Re<Jdc(T) (cf. Fig. |Hl lower panel). Since this 
signature is related to (rather local) polaron excitation 
processes the position of the maximum is almost inde- 
pendent of the system size. In comparison to the d = oo 
(DMFT) results^ we find the same qualitative behav- 
ior in the relevant temperature regime loq < T ~ 2e p 
but a different location of the conductivity maximum. 
Generally in DMFT the activation energy for polaron 
hopping turns out to be lower than expected from com- 
monly accepted arguments for finite-d systems. Increas- 
ing the lattice size, our Id Holstein data indicates that 
this discrepancy does not necessarily imply the failure of 
standard theory of hopping conduction^ but may partly 
arise from dimensionality effects on polaron transport in 
infinite dimensions. Conversely, and in light of the de- 
viations found for the ac conductivity (cf. Fig. 0J, the 
standard (anti-adiabatic) strong-coupling description can 
only be supposed to provide estimates on relevant energy 
scales in the intermediate adiabatic EP-coupling regime. 



„ D-l 

Rea dc = || e ~^ En IHJHI 2 S(E m - E n ) . (11) 

Re a ^ essentially counts the number of thermally acces- 
sible current carrying (degenerate) states. Since we have 
(0|J|0) = the conductivity almost vanishes for small T 
below the finite size gap between the ground state and the 
first excited state. In the thermodynamic limit, L — > oo, 
Re<7dc is related to the charge stiffness D c , or the so- 
called Drude weight (at T = 0)M 

The temperature dependence of the dc conductivity is 
illustrated in Fig. |S| Again the weak coupling results ap- 
pear to have a rather strong finite-size dependence (note 
that Re <7dc depicted in Fig. |S1 is an intensive quantity) . 
When L increases a pronounced peak develops at low 
temperatures. This peak can be attributed to a nor- 
mal "metallic-like" behavior. Independent of the cou- 
pling strength e p the polaron is an itinerant quasiparti- 
cle thus, for T — > 0, always leading to band conduction. 
When e p > 0, we expect that Reer(w — > 0) is finite for 
T > in contrast to an ideal conductor. 

The shoulder observed for the 10 and 12 site systems 
at T/t > 1, again is an artifact of our phonon trunca- 
tion procedure as can be seen by comparing the data 
obtained for L = 6 and different numbers of the phonon 
cut-off M. In the strong-EP-coupling polaronic regime, 
band-like transport becomes extensively suppressed (the 
Drude weight is exponentially small). Nevertheless quan- 
tum zero-point phonon fluctuations cause polaron der- 
ealization at T = 0. At higher temperatures incoherent 



IV. SUMMARY 

In this work, we have investigated the motion of a 
charge carrier in response to ac and dc external fields 
for strongly correlated Id electron-phonon systems. The 
combination of Lanczos diagonalization and the Kernel 
Polynomial Method has enabled us to calculate for the 
first time quasi-exactly the temperature dependence of 
the optical absorption spectra and the dc conductivity 
in the framework of the one-dimensional Holstein model. 
Besides the well-known polaron maximum a pronounced 
absorption feature at about 2t is found in the optical con- 
ductivity. Finite-size effects were identified and assessed, 
e.g., on the basis of the /-sum rule. In the physically most 
interesting range of intermediate coupling strengths and 
phonon frequencies, we find that the conductivity devi- 
ates from the standard small polaron results. 
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